(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 12.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     73752,       1735]
NotebookOptionsPosition[     70264,       1681]
NotebookOutlinePosition[     70650,       1698]
CellTagsIndexPosition[     70607,       1695]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:6e90\:7a0b\:5e8f\:8fc7\:4e8e\:590d\:6742\:548c\:5197\:957f", "\:ff0c", 
    "\:6ca1\:6cd5\:770b", "\:ff0c", "\:91cd\:5199"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"\:76ee\:6807", ":", "\:8f93\:5165\:4e24\:4e2a\:8fb9\:957f"}], 
    "\:ff1a", 
    RowBox[{
    "l3", " ", "l4", " ", 
     "\:8f93\:51fa\:516d\:4e2a\:8282\:70b9\:5750\:6807"}], "\:ff1a", 
    RowBox[{"{", 
     RowBox[{"v1", ",", "v3", ",", "v5", ",", "v2", ",", "v4", ",", "v6"}], 
     "}"}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"\:8981\:6c42", ":", "1"}], ")"}], " ", 
     "\:4e3b\:51fd\:6570Vertices", " ", 
     "\:53ef\:4ee5\:53e6\:5b58\:7a0b\:5e8f\:5305", " ", 
     RowBox[{"vertice", ".", "wdx"}]}], "\:ff0c", 
    "\:53ef\:4ee5\:5355\:72ec\:4e8e\:5176\:4ed6\:7a0b\:5e8f\:8fd0\:884c"}], 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"     ", "2", ")"}], " ", 
    "\:4e3b\:51fd\:6570\:7b26\:5408\:67f1\:578b\:5bf9\:79f0\:6761\:4ef6"}], 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"     ", "3", ")"}], " ", 
    "\:8f93\:51fa\:8282\:70b9\:6784\:6210simplicial", " ", "complex", " ", 
    "\:7b26\:5408\:4e09\:89d2\:4e0d\:7b49\:5f0f\:548c\:7279\:5b9a\:7684\:5b9a\
\:5411"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"--", 
    RowBox[{"--", 
     RowBox[{"--", 
      RowBox[{"--", 
       RowBox[{"--", 
        RowBox[{"--", 
         RowBox[{"--", 
          RowBox[{"--", 
           RowBox[{"--", 
            RowBox[{"--", 
             RowBox[{"--", 
              RowBox[{"--", 
               RowBox[{"--", 
                RowBox[{"--", 
                 RowBox[{"--", 
                  RowBox[{"--", 
                   RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{
                    "--", "--"}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]\
}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"\:8f85\:52a9\:51fd\:65701", "\:ff1a", " ", 
    RowBox[{"\:8ddd\:79bb", " ", 
     RowBox[{"MVdistance", "[", 
      RowBox[{"VR_", ",", "a_", ",", "b_"}], "]"}], "  ", "\:8f93\:5165"}], 
    "\:ff1a", 
    RowBox[{"VR", " ", "\:516d\:8282\:70b9\:5750\:6807"}], "\:ff0c", " ", "a",
     "\:ff0c", 
    RowBox[{"b", " ", "\:8282\:70b9\:7f16\:53f7", " ", "\:8f93\:51fa"}], 
    "\:ff1a", 
    RowBox[{"\:95f5\:6c0f\:8ddd\:79bb", " ", "\:4f9d\:8d56"}], "\:ff1a", 
    "\:65e0"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"\:8f85\:52a9\:51fd\:65702", "\:ff1a", " ", 
    RowBox[{
     RowBox[{"Vol", "[", 
      RowBox[{"VR_", ",", 
       RowBox[{"{", 
        RowBox[{"a_", ",", "b_", ",", "c_", ",", "..."}], "}"}], ",", "E"}], 
      "]"}], "  ", "\:8f93\:5165"}], "\:ff1a", 
    RowBox[{"VR", " ", "\:516d\:8282\:70b9\:5750\:6807"}], "\:ff0c", " ", "a",
     "\:ff0c", "b", ",", 
    RowBox[{
     RowBox[{
     "...", " ", "\:8282\:70b9\:7f16\:53f7", " ", "E", " ", "\:6b27\:5f0f"}], 
     "=", "1"}], "\:ff0c", 
    RowBox[{"\:6d1b\:4f26\:5179", "=", 
     RowBox[{"0", " ", "\:8f93\:51fa"}]}], "\:ff1a", 
    RowBox[{"n\:7ef4\:4f53\:79ef", " ", "\:4f9d\:8d56"}], "\:ff1a", "CMM"}], 
   " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"\:8f85\:52a9\:51fd\:65702", "\:ff1a", " ", 
    RowBox[{
     RowBox[{"CMM", "[", 
      RowBox[{"VR_", ",", 
       RowBox[{"{", 
        RowBox[{"a_", ",", "b_", ",", "c_", ",", "..."}], "}"}]}], "]"}], "  ",
      "\:8f93\:5165"}], "\:ff1a", 
    RowBox[{"VR", " ", "\:516d\:8282\:70b9\:5750\:6807"}], "\:ff0c", " ", "a",
     "\:ff0c", "b", ",", 
    RowBox[{"...", " ", "\:8282\:70b9\:7f16\:53f7", " ", "\:8f93\:51fa"}], 
    "\:ff1a", 
    RowBox[{"n", " ", "Caley", " ", "Matrix", " ", "\:4f9d\:8d56"}], "\:ff1a",
     "\:65e0"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"--", 
    RowBox[{"--", 
     RowBox[{"--", 
      RowBox[{"--", 
       RowBox[{"--", 
        RowBox[{"--", 
         RowBox[{"--", 
          RowBox[{"--", 
           RowBox[{"--", 
            RowBox[{"--", 
             RowBox[{"--", 
              RowBox[{"--", 
               RowBox[{"--", 
                RowBox[{"--", 
                 RowBox[{"--", 
                  RowBox[{"--", 
                   RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{
                    "--", "--"}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]\
}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"\:4e3b\:8981\:6570\:636e\:63a5\:53e3", ":", " ", "1"}], ")"}], 
    " ", "VR", " ", "\:516d\:8282\:70b9\:5750\:6807", " ", 
    RowBox[{"(", "v6\:672a\:786e\:5b9a\:524d\:4e3a0", ")"}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"               ", "2", ")"}], " ", "vid", "  ", "VR", " ", 
    "\:4e2d\:7ed3\:70b9\:7f16\:53f7", " ", 
    RowBox[{"(", "v6\:672a\:786e\:5b9a\:524d\:4e3a0", ")"}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"               ", "3", ")"}], " ", 
     RowBox[{"SYM", "[", 
      RowBox[{"[", 
       RowBox[{"inputs", ",", "edges", ",", "vertices"}], "]"}], "]"}], " ", 
     "l1"}], ",", "l2", ",", "l3", ",", "l4", ",", 
    "l5\:5bf9\:5e94\:7684edge"}], "*)"}], "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.8240330633124676`*^9, 3.8240331046389475`*^9}, {
  3.8240760608112636`*^9, 3.8240760637007256`*^9}, {3.8240761113780165`*^9, 
  3.82407616691424*^9}, {3.824076213330163*^9, 3.8240765652215614`*^9}, {
  3.824076599169849*^9, 3.8240766677220197`*^9}, {3.8240768421432824`*^9, 
  3.824076843951989*^9}, {3.8240786038367953`*^9, 3.824078622270731*^9}, {
  3.8240786655711813`*^9, 3.8240788079976416`*^9}, {3.82407902758671*^9, 
  3.8240791140449457`*^9}, {3.8240793922845144`*^9, 3.8240795348927245`*^9}, {
  3.8240942549244757`*^9, 3.8240943841245704`*^9}, {3.8240949925116296`*^9, 
  3.8240949969186883`*^9}, {3.8241058881729937`*^9, 3.8241058943466487`*^9}, {
  3.824106106353774*^9, 3.8241061999885283`*^9}, {3.824158551568694*^9, 
  3.824158555865343*^9}},ExpressionUUID->"398a21c6-c9a3-4ee1-8c7d-\
444260a91e72"],

Cell[BoxData[
 RowBox[{"Clear", "[", "\"\<Global`*\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.817496800058392*^9, 3.817496801305118*^9}, 
   3.824012658387802*^9},
 CellLabel->
  "In[355]:=",ExpressionUUID->"71db4609-c829-4834-b1c4-1c1feed1af1c"],

Cell[BoxData[
 RowBox[{
  RowBox[{"SetDirectory", "[", 
   RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.824012670220211*^9, 3.8240126766254797`*^9}, {
  3.8240130924222403`*^9, 3.824013093975934*^9}},
 CellLabel->
  "In[357]:=",ExpressionUUID->"e1070aa4-d539-4b8c-b081-7402b6b16374"],

Cell[BoxData[
 RowBox[{
  RowBox[{"tempfile13", "=", 
   RowBox[{"FileNameJoin", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Directory", "[", "]"}], ",", "\"\<vertice.wdx\>\""}], "}"}], 
    "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.8240126898708444`*^9, 3.824012691646351*^9}, {
  3.824013028331542*^9, 3.8240130299380026`*^9}},
 CellLabel->
  "In[358]:=",ExpressionUUID->"d70282d0-7995-4523-b08c-3a45a4e5f7a1"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"vid", "=", 
    RowBox[{"{", 
     RowBox[{"1", ",", "3", ",", "5", ",", "2", ",", "4", ",", "6"}], "}"}]}],
    ";"}], 
  RowBox[{"(*", 
   RowBox[{"order", " ", "of", " ", "vertices"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.8240367179852905`*^9, 3.8240367296134872`*^9}, {
  3.8240768535038395`*^9, 3.824076871245906*^9}},
 CellLabel->
  "In[359]:=",ExpressionUUID->"58628041-6b46-49e7-b651-4abb3f7a373a"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"SYM", "=", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"1", ",", "2"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"1", ",", "3"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"2", ",", "3"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"1", ",", "5"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"3", ",", "5"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"2", ",", "5"}], "}"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"3", ",", "4"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"1", ",", "4"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"2", ",", "4"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"4", ",", "5"}], "}"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"1", ",", "6"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"3", ",", "6"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"5", ",", "6"}], "}"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"{", 
        RowBox[{"2", ",", "6"}], "}"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"{", 
        RowBox[{"4", ",", "6"}], "}"}], "}"}]}], "}"}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"[", 
    RowBox[{"[", 
     RowBox[{"inputs", ",", "edges", ",", "vertices"}], "]"}], "]"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Sid", "=", 
   RowBox[{"{", 
    RowBox[{"2", ",", "4", ",", "6"}], "}"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.8240768111625214`*^9, 3.824076830308875*^9}, {
  3.824076895049386*^9, 3.8240769055459967`*^9}, {3.824104827692011*^9, 
  3.8241048421226134`*^9}},
 CellLabel->
  "In[360]:=",ExpressionUUID->"e80a3ff6-3336-4d4f-a00d-95c9256a2780"],

Cell[BoxData[
 RowBox[{
  RowBox[{"ge", "=", 
   RowBox[{"DiagonalMatrix", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "1"}], ",", "1", ",", "1", ",", "1"}], "}"}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.817483534499908*^9, 3.817483539645171*^9}, {
   3.824035629106094*^9, 3.8240356304066906`*^9}, 3.824103503895243*^9, 
   3.824103672393691*^9},
 CellLabel->
  "In[362]:=",ExpressionUUID->"4c998990-0f46-443b-ac3a-a12e4bd34016"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Vertex", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", 
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"2", " ", 
            RowBox[{"\[Sqrt]", "5"}]}], ")"}], "/", 
          SuperscriptBox["3", 
           RowBox[{"1", "/", "4"}]]}], ")"}]}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", 
       RowBox[{
        RowBox[{"-", 
         SuperscriptBox["3", 
          RowBox[{"1", "/", "4"}]]}], " ", 
        RowBox[{"\[Sqrt]", "5"}]}], ",", 
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"\[Sqrt]", "5"}], ")"}], "/", 
          SuperscriptBox["3", 
           RowBox[{"1", "/", "4"}]]}], ")"}]}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"2", " ", 
            RowBox[{"\[Sqrt]", "10"}]}], ")"}], "/", 
          SuperscriptBox["3", 
           RowBox[{"3", "/", "4"}]]}], ")"}]}], ",", 
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"\[Sqrt]", "5"}], ")"}], "/", 
          SuperscriptBox["3", 
           RowBox[{"3", "/", "4"}]]}], ")"}]}], ",", 
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"\[Sqrt]", "5"}], ")"}], "/", 
          SuperscriptBox["3", 
           RowBox[{"1", "/", "4"}]]}], ")"}]}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{"1", "/", 
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox["3", 
             RowBox[{"1", "/", "4"}]], " ", 
            RowBox[{"\[Sqrt]", "10"}]}], ")"}]}], ")"}]}], ",", 
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"\[Sqrt]", 
            RowBox[{"(", 
             RowBox[{"5", "/", "2"}], ")"}]}], ")"}], "/", 
          SuperscriptBox["3", 
           RowBox[{"3", "/", "4"}]]}], ")"}]}], ",", 
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"\[Sqrt]", "5"}], ")"}], "/", 
          SuperscriptBox["3", 
           RowBox[{"3", "/", "4"}]]}], ")"}]}], ",", 
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"\[Sqrt]", "5"}], ")"}], "/", 
          SuperscriptBox["3", 
           RowBox[{"1", "/", "4"}]]}], ")"}]}]}], "}"}]}], "}"}]}], ";", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"vertices", " ", "v1"}], "\[Rule]", " ", 
    RowBox[{"v5", " ", "from", " ", "previous", " ", "computation"}]}], 
   "*)"}]}]], "Input",
 CellChangeTimes->{{3.824077919989416*^9, 3.8240779569325943`*^9}, 
   3.8241070499760647`*^9},
 CellLabel->
  "In[363]:=",ExpressionUUID->"68b8b8fb-cf6f-4209-98ef-b208f8ed33c4"],

Cell[BoxData[
 RowBox[{
  RowBox[{"MVdistance", "[", 
   RowBox[{"VR_", ",", "a_", ",", "b_"}], "]"}], ":=", 
  RowBox[{"Block", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"ge", "=", 
      RowBox[{"DiagonalMatrix", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "1"}], ",", "1", ",", "1", ",", "1"}], "}"}], "]"}]}], 
     "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"VR", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"Position", "[", 
           RowBox[{"vid", ",", "a"}], "]"}], "[", 
          RowBox[{"[", 
           RowBox[{"1", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", 
       RowBox[{"VR", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"Position", "[", 
           RowBox[{"vid", ",", "b"}], "]"}], "[", 
          RowBox[{"[", 
           RowBox[{"1", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], ".", 
     "ge", ".", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"VR", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"Position", "[", 
           RowBox[{"vid", ",", "a"}], "]"}], "[", 
          RowBox[{"[", 
           RowBox[{"1", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", 
       RowBox[{"VR", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"Position", "[", 
           RowBox[{"vid", ",", "b"}], "]"}], "[", 
          RowBox[{"[", 
           RowBox[{"1", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}]}], 
   "\[IndentingNewLine]", "]"}]}]], "Input",
 CellLabel->
  "In[364]:=",ExpressionUUID->"d2601902-6df0-4653-a0ef-ed3cfe9da0ae"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"SOL", "[", 
    RowBox[{"Vertex_", ",", "l3_", ",", "l4_"}], "]"}], ":=", 
   RowBox[{"Block", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"VK", ",", "Sol", ",", "t", ",", "x", ",", "y", ",", "z"}], 
      "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"VK", "=", 
       RowBox[{"Append", "[", 
        RowBox[{"Vertex", ",", 
         RowBox[{"{", 
          RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Sol", "=", 
       RowBox[{"Solve", "[", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"MVdistance", "[", 
            RowBox[{"VK", ",", 
             RowBox[{"SYM", "[", 
              RowBox[{"[", 
               RowBox[{"4", ",", "1", ",", "1"}], "]"}], "]"}], ",", 
             RowBox[{"SYM", "[", 
              RowBox[{"[", 
               RowBox[{"4", ",", "1", ",", "2"}], "]"}], "]"}]}], "]"}], 
           "\[Equal]", 
           RowBox[{"l4", "^", "2"}]}], "&&", 
          RowBox[{
           RowBox[{"MVdistance", "[", 
            RowBox[{"VK", ",", 
             RowBox[{"SYM", "[", 
              RowBox[{"[", 
               RowBox[{"3", ",", "1", ",", "1"}], "]"}], "]"}], ",", 
             RowBox[{"SYM", "[", 
              RowBox[{"[", 
               RowBox[{"3", ",", "1", ",", "2"}], "]"}], "]"}]}], "]"}], 
           "\[Equal]", 
           RowBox[{"l3", "^", "2"}]}], "&&", 
          RowBox[{
           RowBox[{"MVdistance", "[", 
            RowBox[{"VK", ",", 
             RowBox[{"SYM", "[", 
              RowBox[{"[", 
               RowBox[{"3", ",", "2", ",", "1"}], "]"}], "]"}], ",", 
             RowBox[{"SYM", "[", 
              RowBox[{"[", 
               RowBox[{"3", ",", "2", ",", "2"}], "]"}], "]"}]}], "]"}], 
           "\[Equal]", 
           RowBox[{"l3", "^", "2"}]}], "&&", 
          RowBox[{
           RowBox[{"MVdistance", "[", 
            RowBox[{"VK", ",", 
             RowBox[{"SYM", "[", 
              RowBox[{"[", 
               RowBox[{"3", ",", "3", ",", "1"}], "]"}], "]"}], ",", 
             RowBox[{"SYM", "[", 
              RowBox[{"[", 
               RowBox[{"3", ",", "3", ",", "2"}], "]"}], "]"}]}], "]"}], 
           "\[Equal]", 
           RowBox[{"l3", "^", "2"}]}]}], ",", 
         RowBox[{"{", 
          RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Map", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"VK", "/.", "#"}], "&"}], ",", "Sol"}], "]"}]}]}], 
    "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.8240958754571447`*^9, 3.8240959035554795`*^9}, {
  3.8241026361341977`*^9, 3.8241027910507116`*^9}, {3.8241072484247894`*^9, 
  3.824107306208644*^9}, {3.824107413763712*^9, 3.824107450919898*^9}},
 CellLabel->
  "In[365]:=",ExpressionUUID->"ad538921-e09d-4dd7-a1d4-c57ba0f92e3f"],

Cell[BoxData[
 RowBox[{
  RowBox[{"CMM", "[", 
   RowBox[{"VK_", ",", "x_"}], "]"}], ":=", 
  RowBox[{"Block", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"N", "=", 
       RowBox[{"Length", "[", "x", "]"}]}], ",", "DD"}], "}"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"DD", "=", 
      RowBox[{"PadLeft", "[", 
       RowBox[{
        RowBox[{"Table", "[", 
         RowBox[{
          RowBox[{"MVdistance", "[", 
           RowBox[{"VK", ",", 
            RowBox[{"x", "[", 
             RowBox[{"[", "i", "]"}], "]"}], ",", 
            RowBox[{"x", "[", 
             RowBox[{"[", "j", "]"}], "]"}]}], "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"i", ",", "1", ",", "N"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"j", ",", "1", ",", "N"}], "}"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"N", "+", "1"}], ",", 
          RowBox[{"N", "+", "1"}]}], "}"}], ",", "1"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"DD", "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "1"}], "]"}], "]"}], "=", "0"}], ";", 
     "\[IndentingNewLine]", "DD"}]}], "\[IndentingNewLine]", "]"}]}]], "Input",\

 CellChangeTimes->{{3.8241026194638987`*^9, 3.8241026311047697`*^9}, {
   3.8241028861446824`*^9, 3.824103042260583*^9}, 3.824103099980062*^9, {
   3.8241036384806023`*^9, 3.824103641159957*^9}},
 CellLabel->
  "In[366]:=",ExpressionUUID->"19b11aa4-8a24-438f-b35a-d67d7ad30830"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Vol", "[", 
   RowBox[{"VK_", ",", "x_", ",", "E_"}], "]"}], ":=", 
  RowBox[{"Block", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"DD", ",", 
      RowBox[{"N", "=", 
       RowBox[{
        RowBox[{"Length", "[", "x", "]"}], "-", "1"}]}]}], "}"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"DD", "=", 
      RowBox[{"CMM", "[", 
       RowBox[{"VK", ",", "x"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"Sqrt", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"-", "1"}], ")"}], "^", 
         RowBox[{"(", 
          RowBox[{"N", "+", "E"}], ")"}]}], "/", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"2", "^", "N"}], "*", 
          RowBox[{
           RowBox[{"Factorial", "[", "N", "]"}], "^", "2"}]}], ")"}]}], "*", 
       RowBox[{"Det", "[", "DD", "]"}]}], "]"}]}]}], "\[IndentingNewLine]", 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.8241058686015863`*^9, 3.824106007780486*^9}, {
  3.8241060675060606`*^9, 3.8241060688828697`*^9}},
 CellLabel->
  "In[367]:=",ExpressionUUID->"d0e2b5ef-98e0-4c6f-81b9-ad4ff82e8a0d"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Vertices", "[", 
   RowBox[{"vid_", ",", "Sid_", ",", "Vertex_", ",", "l3_", ",", "l4_"}], 
   "]"}], ":=", 
  RowBox[{"Block", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"VK", ",", "A"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"A", "=", 
      RowBox[{"Map", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"vid", "/.", 
          RowBox[{"#", "\[Rule]", "Nothing"}]}], "&"}], ",", "Sid"}], "]"}]}],
      ";", "\[IndentingNewLine]", 
     RowBox[{"VK", "=", 
      RowBox[{
       RowBox[{"SOL", "[", 
        RowBox[{"Vertex", ",", "l3", ",", "l4"}], "]"}], "//", "N"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{
          RowBox[{"Element", "[", 
           RowBox[{
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"Map", "[", 
               RowBox[{
                RowBox[{
                 RowBox[{"Vol", "[", 
                  RowBox[{
                   RowBox[{"VK", "[", 
                    RowBox[{"[", "z", "]"}], "]"}], ",", "#", ",", "1"}], 
                  "]"}], "&"}], ",", 
                RowBox[{"Subsets", "[", 
                 RowBox[{
                  RowBox[{"A", "[", 
                   RowBox[{"[", "x", "]"}], "]"}], ",", 
                  RowBox[{"{", "3", "}"}]}], "]"}]}], "]"}], ",", 
              RowBox[{"{", 
               RowBox[{"x", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", 
            "Reals"}], "]"}], "&&", 
          RowBox[{"Element", "[", 
           RowBox[{
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"Map", "[", 
               RowBox[{
                RowBox[{
                 RowBox[{"Vol", "[", 
                  RowBox[{
                   RowBox[{"VK", "[", 
                    RowBox[{"[", "z", "]"}], "]"}], ",", "#", ",", "1"}], 
                  "]"}], "&"}], ",", 
                RowBox[{"Subsets", "[", 
                 RowBox[{
                  RowBox[{"A", "[", 
                   RowBox[{"[", "x", "]"}], "]"}], ",", 
                  RowBox[{"{", "4", "}"}]}], "]"}]}], "]"}], ",", 
              RowBox[{"{", 
               RowBox[{"x", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", 
            "Reals"}], "]"}], "&&", 
          RowBox[{"Element", "[", 
           RowBox[{
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"Map", "[", 
               RowBox[{
                RowBox[{
                 RowBox[{"Vol", "[", 
                  RowBox[{
                   RowBox[{"VK", "[", 
                    RowBox[{"[", "z", "]"}], "]"}], ",", "#", ",", "0"}], 
                  "]"}], "&"}], ",", 
                RowBox[{"Subsets", "[", 
                 RowBox[{
                  RowBox[{"A", "[", 
                   RowBox[{"[", "x", "]"}], "]"}], ",", 
                  RowBox[{"{", "5", "}"}]}], "]"}]}], "]"}], ",", 
              RowBox[{"{", 
               RowBox[{"x", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", 
            "Reals"}], "]"}]}], ",", 
         RowBox[{"VK", "[", 
          RowBox[{"[", "z", "]"}], "]"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"z", ",", "1", ",", 
         RowBox[{"Length", "[", "VK", "]"}]}], "}"}]}], "]"}]}]}], 
   "\[IndentingNewLine]", "]"}]}]], "Input",
 CellChangeTimes->{{3.8241066883586717`*^9, 3.8241068646631885`*^9}, {
  3.8241071645060396`*^9, 3.8241071679280243`*^9}, {3.8241076514763374`*^9, 
  3.824107661194109*^9}},
 CellLabel->
  "In[368]:=",ExpressionUUID->"69850235-1da2-4e23-befc-4c70c9199bbb"],

Cell[BoxData[
 RowBox[{
  RowBox[{"DumpSave", "[", 
   RowBox[{"tempfile13", ",", 
    RowBox[{"{", 
     RowBox[{
     "Vertices", ",", "MVdistance", ",", "Vol", ",", "CMM", ",", "vid", ",", 
      "Vertex", ",", "Sid", ",", "SOL", ",", "SYM"}], "}"}]}], "]"}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.824013043775469*^9, 3.8240130619554377`*^9}, {
  3.8240139677884846`*^9, 3.8240139755411873`*^9}, {3.8241078642930984`*^9, 
  3.824107917022148*^9}, {3.824108141219358*^9, 3.824108143559636*^9}, {
  3.824108192223736*^9, 3.824108193296647*^9}},
 CellLabel->
  "In[369]:=",ExpressionUUID->"6c01a0ed-c53e-4a04-ac7b-64a4bf5a41d4"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Vertices", "[", 
   RowBox[{"vid", ",", "Sid", ",", "Vertex", ",", 
    RowBox[{"Sqrt", "[", "10", "]"}], ",", 
    RowBox[{"Sqrt", "[", "31", "]"}]}], "]"}], "[", 
  RowBox[{"[", "1", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.8241068700326357`*^9, 3.824106899095076*^9}, {
  3.8241070974871264`*^9, 3.8241071122469997`*^9}, {3.824107627755985*^9, 
  3.824107674131596*^9}, {3.824107784184263*^9, 3.8241078345857377`*^9}},
 CellLabel->
  "In[370]:=",ExpressionUUID->"d447adda-a031-4fc8-8a5b-5fc66c8cc0b8"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"True", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"0.`", ",", "0.`", ",", "0.`", ",", "0.`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.`", ",", "0.`", ",", "0.`", ",", 
       RowBox[{"-", "3.398088489694245`"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.`", ",", "0.`", ",", 
       RowBox[{"-", "2.942830956382712`"}], ",", 
       RowBox[{"-", "1.6990442448471226`"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.`", ",", 
       RowBox[{"-", "2.7745276335252114`"}], ",", 
       RowBox[{"-", "0.9809436521275706`"}], ",", 
       RowBox[{"-", "1.6990442448471226`"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "0.24028114141347542`"}], ",", 
       RowBox[{"-", "0.6936319083813028`"}], ",", 
       RowBox[{"-", "0.9809436521275706`"}], ",", 
       RowBox[{"-", "1.6990442448471226`"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "1.8444571380239085`"}], ",", "3.0907960688809353`", ",", 
       
       RowBox[{"-", "0.9809436521275706`"}], ",", 
       RowBox[{"-", "1.6990442448471226`"}]}], "}"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.8241068998689795`*^9, 3.8241070624621067`*^9, {3.8241071006130333`*^9, 
   3.8241071134567127`*^9}, 3.824107171344146*^9, {3.8241076320339537`*^9, 
   3.8241076748986597`*^9}, {3.8241077850939903`*^9, 3.824107835507191*^9}, 
   3.824107947870388*^9, 3.8241081693593674`*^9, 3.824108226099625*^9},
 CellLabel->
  "Out[370]=",ExpressionUUID->"4174b2cb-305b-4658-bc33-7c5eec2aeeba"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"Areas", "[", 
   RowBox[{"l1_", ",", "l2_", ",", "l3_"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"s", ",", "tarea"}], "}"}], ",", " ", 
    RowBox[{
     RowBox[{"s", "=", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"l1", "+", "l2", "+", "l3"}], ")"}], "/", "2"}]}], ";", 
     RowBox[{"tarea", "=", 
      RowBox[{"Sqrt", "[", 
       RowBox[{"s", 
        RowBox[{"(", 
         RowBox[{"s", "-", "l1"}], ")"}], 
        RowBox[{"(", 
         RowBox[{"s", "-", "l2"}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{"s", "-", "l3"}], ")"}]}], "]"}]}], ";", " ", 
     RowBox[{"Return", "[", "tarea", "]"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.823889460177886*^9, 3.823889489137599*^9}, {
  3.823889541215632*^9, 3.8238895644530573`*^9}, {3.823889612977397*^9, 
  3.8238896775952063`*^9}, {3.823889724330112*^9, 3.8238897261383963`*^9}, {
  3.82388978491133*^9, 3.823889789965691*^9}},
 CellLabel->"In[58]:=",ExpressionUUID->"4f3c4822-5dd2-4a93-9e50-e5e37aa00ae9"],

Cell[BoxData[
 RowBox[{
  RowBox[{"AreaTriangle", "[", 
   RowBox[{"{", 
    RowBox[{"x1_", ",", "x2_", ",", "x3_"}], "}"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"l1", ",", "l2", ",", "tarea"}], "}"}], ",", " ", 
    RowBox[{
     RowBox[{"l1", "=", 
      RowBox[{"x1", "-", "x2"}]}], ";", " ", 
     RowBox[{"l2", "=", 
      RowBox[{"x1", "-", "x3"}]}], ";", " ", 
     RowBox[{"tarea", "=", 
      RowBox[{
       RowBox[{"1", "/", "2"}], 
       RowBox[{"Sqrt", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{"l1", ".", "\[Eta]", ".", "l1"}], ")"}], 
          RowBox[{"(", 
           RowBox[{"l2", ".", "\[Eta]", ".", "l2"}], ")"}]}], " ", "-", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"l1", ".", "\[Eta]", ".", "l2"}], ")"}], "^", "2"}]}], 
        "]"}]}]}], ";", " ", 
     RowBox[{"Return", "[", "tarea", "]"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.696770053562324*^9, 3.696770053824493*^9}, {
   3.6968490227161016`*^9, 3.6968490460799437`*^9}, {3.696849099775785*^9, 
   3.696849125512746*^9}, {3.69684923914423*^9, 3.696849247727893*^9}, {
   3.699162534214698*^9, 3.699162539049659*^9}, {3.6991625869453993`*^9, 
   3.699162615377471*^9}, 3.699190648688534*^9, {3.699248357103242*^9, 
   3.69924835758105*^9}, {3.6992556838335943`*^9, 3.699255697080563*^9}, {
   3.699255963028359*^9, 3.699255968160611*^9}, {3.699269447030582*^9, 
   3.699269457054044*^9}, {3.699939266867931*^9, 3.699939276209908*^9}, {
   3.699939907040021*^9, 3.6999399114350033`*^9}, {3.6999420267079573`*^9, 
   3.699942029263701*^9}, {3.73412251529758*^9, 3.7341225522742367`*^9}, {
   3.734122605823371*^9, 3.734122608215273*^9}, {3.7341250444245167`*^9, 
   3.734125067024579*^9}, {3.734125100853516*^9, 3.7341251263756027`*^9}, {
   3.7341251699894037`*^9, 3.734125186987216*^9}, {3.734125876631937*^9, 
   3.7341259114254*^9}, {3.734128174601564*^9, 3.7341282089722977`*^9}, {
   3.7341773716009083`*^9, 3.734177422411697*^9}, {3.73418600765477*^9, 
   3.734186019550564*^9}, {3.734449615409775*^9, 3.734449676434636*^9}, {
   3.734449771521804*^9, 3.7344497769860153`*^9}, {3.734450093262422*^9, 
   3.7344500951356783`*^9}, {3.748197742304473*^9, 3.748197979019801*^9}, {
   3.7481980631403017`*^9, 3.748198104518581*^9}, {3.748198192134512*^9, 
   3.748198194421586*^9}, 3.748198497595421*^9, {3.748198717586803*^9, 
   3.748198794413993*^9}, {3.748200742389719*^9, 3.7482007484437113`*^9}, {
   3.748201365965706*^9, 3.7482013664993753`*^9}, {3.817573815958921*^9, 
   3.817573816632172*^9}},
 CellLabel->
  "In[108]:=",ExpressionUUID->"64b89cc3-7eae-45f0-9067-6ee016e76151"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Tetnormal", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"x1_", ",", "x2_", ",", "x3_", ",", "x4_"}], "}"}], ",", "c_"}], 
   "]"}], ":=", " ", 
  RowBox[{"Block", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"l1", ",", "l2", ",", "l3", ",", "nn", ",", "norm"}], "}"}], ",", 
    RowBox[{
     RowBox[{"l1", "=", 
      RowBox[{"x1", "-", "x2"}]}], ";", " ", 
     RowBox[{"l2", "=", 
      RowBox[{"x1", "-", "x3"}]}], ";", 
     RowBox[{"l3", "=", 
      RowBox[{"x1", "-", "x4"}]}], ";", " ", 
     RowBox[{"nn", "=", 
      RowBox[{
       RowBox[{"LeviCivitaTensor", "[", "4", "]"}], ".", 
       RowBox[{"(", 
        RowBox[{"\[Eta]", ".", "l1"}], ")"}], ".", 
       RowBox[{"(", 
        RowBox[{"\[Eta]", ".", "l2"}], ")"}], ".", 
       RowBox[{"(", 
        RowBox[{"\[Eta]", ".", "l3"}], ")"}]}]}], ";", " ", 
     RowBox[{"norm", " ", "=", " ", 
      RowBox[{"Sqrt", "[", 
       RowBox[{"-", 
        RowBox[{"nn", ".", "\[Eta]", ".", "nn"}]}], "]"}]}], ";", " ", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"Refine", "[", 
        RowBox[{
         RowBox[{"nn", ".", 
          RowBox[{"(", 
           RowBox[{"\[Eta]", ".", 
            RowBox[{"(", 
             RowBox[{"x1", "-", "c"}], ")"}]}], ")"}]}], ">", "0"}], "]"}], 
       ",", 
       RowBox[{"Return", "[", 
        RowBox[{"nn", "/", "norm"}], "]"}], ",", 
       RowBox[{"Return", "[", 
        RowBox[{
         RowBox[{"-", "nn"}], "/", "norm"}], "]"}]}], "]"}]}]}], " ", 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.818005916228182*^9, 3.818005923663908*^9}, 
   3.818005982009305*^9, {3.818873597719404*^9, 3.8188735984042397`*^9}},
 CellLabel->
  "In[109]:=",ExpressionUUID->"8be7cf6b-277d-4984-9bf5-23ee2193efa9"],

Cell[BoxData[
 RowBox[{
  RowBox[{"CenterT", "[", "tet_", "]"}], ":=", " ", 
  RowBox[{
   RowBox[{"Total", "[", "tet", "]"}], "/", 
   RowBox[{"Length", "[", "tet", "]"}]}]}]], "Input",
 CellLabel->
  "In[110]:=",ExpressionUUID->"890cd635-1b64-4a81-bc66-96214fa89eba"],

Cell[BoxData[
 RowBox[{
  RowBox[{"dihedral", "[", 
   RowBox[{"a_", ",", "b_"}], "]"}], ":=", 
  RowBox[{"If", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"a", ".", "\[Eta]", ".", "b"}], ">", "0"}], ",", 
    RowBox[{"ArcCosh", "[", 
     RowBox[{"a", ".", "\[Eta]", ".", "b"}], "]"}], ",", 
    RowBox[{"ArcCosh", "[", 
     RowBox[{"-", 
      RowBox[{"a", ".", "\[Eta]", ".", "b"}]}], "]"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.798543768743003*^9, 3.79854378118576*^9}, {
   3.8175766941832943`*^9, 3.8175767162995872`*^9}, {3.817576777703725*^9, 
   3.817576834310954*^9}, {3.818007815589532*^9, 3.8180078189505463`*^9}, {
   3.818008421188149*^9, 3.81800842567491*^9}, {3.818008456482274*^9, 
   3.818008459300948*^9}, {3.818009685261856*^9, 3.8180096882168283`*^9}, {
   3.818012139618486*^9, 3.81801214286366*^9}, {3.818012530634139*^9, 
   3.818012533250985*^9}, {3.818013234298192*^9, 3.818013237434477*^9}, 
   3.8180206985622797`*^9, 3.818020845146124*^9, {3.8180211530373898`*^9, 
   3.818021159794265*^9}, {3.818110693563748*^9, 3.818110696632753*^9}, {
   3.818881172322996*^9, 3.8188812294383593`*^9}, 3.818881262132627*^9},
 CellLabel->
  "In[111]:=",ExpressionUUID->"d96a3542-e0b4-4569-9ed6-fb27bc6953a4"],

Cell[BoxData[
 RowBox[{
  RowBox[{"dificitAngle", "[", 
   RowBox[{"random3_", ",", "random4_"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "vertice", ",", "Simplices", ",", "FScenter", ",", "tet1", ",", "tet2", 
      ",", "tet3"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"vertice", "=", 
      RowBox[{"Vertices", "[", 
       RowBox[{"random3", ",", "random4"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"vertice", "[", 
        RowBox[{"[", "1", "]"}], "]"}], ",", 
       RowBox[{
        RowBox[{"Simplices", "=", 
         RowBox[{
          RowBox[{
           RowBox[{
            RowBox[{
            "vertice", "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}],
             "[", 
            RowBox[{"[", "#", "]"}], "]"}], "&"}], "/@", "simplices"}]}], ";", 
        RowBox[{"FScenter", "=", 
         RowBox[{"Table", "[", 
          RowBox[{
           RowBox[{"FullSimplify", "@", 
            RowBox[{"CenterT", "@", 
             RowBox[{
             "Simplices", "\[LeftDoubleBracket]", "i", 
              "\[RightDoubleBracket]"}]}]}], ",", 
           RowBox[{"{", 
            RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"tet1", "=", 
         RowBox[{
          RowBox[{
           RowBox[{"Tetnormal", "[", 
            RowBox[{
             RowBox[{
              RowBox[{
              "vertice", "\[LeftDoubleBracket]", "2", 
               "\[RightDoubleBracket]"}], "[", 
              RowBox[{"[", "#", "]"}], "]"}], ",", 
             RowBox[{
             "FScenter", "\[LeftDoubleBracket]", "1", 
              "\[RightDoubleBracket]"}]}], "]"}], "&"}], "/@", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"1", ",", "2", ",", "3", ",", "5"}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{"1", ",", "3", ",", "4", ",", "5"}], "}"}]}], "}"}]}]}], 
        ";", "\[IndentingNewLine]", 
        RowBox[{"tet2", "=", 
         RowBox[{
          RowBox[{
           RowBox[{"Tetnormal", "[", 
            RowBox[{
             RowBox[{
              RowBox[{
              "vertice", "\[LeftDoubleBracket]", "2", 
               "\[RightDoubleBracket]"}], "[", 
              RowBox[{"[", "#", "]"}], "]"}], ",", 
             RowBox[{
             "FScenter", "\[LeftDoubleBracket]", "2", 
              "\[RightDoubleBracket]"}]}], "]"}], "&"}], "/@", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"1", ",", "2", ",", "3", ",", "5"}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{"1", ",", "3", ",", "5", ",", "6"}], "}"}]}], "}"}]}]}], 
        ";", "\[IndentingNewLine]", 
        RowBox[{"tet3", "=", 
         RowBox[{
          RowBox[{
           RowBox[{"Tetnormal", "[", 
            RowBox[{
             RowBox[{
              RowBox[{
              "vertice", "\[LeftDoubleBracket]", "2", 
               "\[RightDoubleBracket]"}], "[", 
              RowBox[{"[", "#", "]"}], "]"}], ",", 
             RowBox[{
             "FScenter", "\[LeftDoubleBracket]", "3", 
              "\[RightDoubleBracket]"}]}], "]"}], "&"}], "/@", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"1", ",", "3", ",", "4", ",", "5"}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{"1", ",", "3", ",", "5", ",", "6"}], "}"}]}], "}"}]}]}], 
        ";", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"tet1", "\[LeftDoubleBracket]", 
              RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}], "*", 
             RowBox[{"tet1", "\[LeftDoubleBracket]", 
              RowBox[{"2", ",", "1"}], "\[RightDoubleBracket]"}]}], ">", 
            "0"}], ",", 
           RowBox[{"-", 
            RowBox[{"dihedral", "[", 
             RowBox[{
              RowBox[{
              "tet1", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], 
              ",", 
              RowBox[{
              "tet1", "\[LeftDoubleBracket]", "2", 
               "\[RightDoubleBracket]"}]}], "]"}]}], ",", 
           RowBox[{"dihedral", "[", 
            RowBox[{
             RowBox[{
             "tet1", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], 
             ",", 
             RowBox[{
             "tet1", "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}]}],
             "]"}]}], "]"}], "+", 
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"tet2", "\[LeftDoubleBracket]", 
              RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}], "*", 
             RowBox[{"tet2", "\[LeftDoubleBracket]", 
              RowBox[{"2", ",", "1"}], "\[RightDoubleBracket]"}]}], ">", 
            "0"}], ",", 
           RowBox[{"-", 
            RowBox[{"dihedral", "[", 
             RowBox[{
              RowBox[{
              "tet2", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], 
              ",", 
              RowBox[{
              "tet2", "\[LeftDoubleBracket]", "2", 
               "\[RightDoubleBracket]"}]}], "]"}]}], ",", 
           RowBox[{"dihedral", "[", 
            RowBox[{
             RowBox[{
             "tet2", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], 
             ",", 
             RowBox[{
             "tet2", "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}]}],
             "]"}]}], "]"}], "+", 
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"tet3", "\[LeftDoubleBracket]", 
              RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}], "*", 
             RowBox[{"tet3", "\[LeftDoubleBracket]", 
              RowBox[{"2", ",", "1"}], "\[RightDoubleBracket]"}]}], ">", 
            "0"}], ",", 
           RowBox[{"-", 
            RowBox[{"dihedral", "[", 
             RowBox[{
              RowBox[{
              "tet3", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], 
              ",", 
              RowBox[{
              "tet3", "\[LeftDoubleBracket]", "2", 
               "\[RightDoubleBracket]"}]}], "]"}]}], ",", 
           RowBox[{"dihedral", "[", 
            RowBox[{
             RowBox[{
             "tet3", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], 
             ",", 
             RowBox[{
             "tet3", "\[LeftDoubleBracket]", "2", "\[RightDoubleBracket]"}]}],
             "]"}]}], "]"}]}]}], ",", "\[IndentingNewLine]", 
       RowBox[{"Continue", "[", "]"}]}], "\[IndentingNewLine]", "]"}]}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.81887504813857*^9, 3.818875080322406*^9}, 
   3.818875897461025*^9, {3.818876001269291*^9, 3.8188760052496557`*^9}, {
   3.81887632046833*^9, 3.818876424967783*^9}, {3.8188767535002947`*^9, 
   3.8188767727266483`*^9}, {3.818876944619138*^9, 3.818876944750209*^9}, {
   3.8188786897337437`*^9, 3.818878748475583*^9}, {3.8188788510521717`*^9, 
   3.8188789357308903`*^9}, {3.818879033359232*^9, 3.818879074579607*^9}, {
   3.818879465291905*^9, 3.818879472300683*^9}, {3.818879524604616*^9, 
   3.8188796585895443`*^9}, {3.818879695252246*^9, 3.8188797160676327`*^9}, 
   3.818879845474187*^9, {3.818880221147793*^9, 3.818880233599408*^9}, {
   3.818880296444994*^9, 3.818880430780451*^9}, {3.818880470658578*^9, 
   3.818880477431913*^9}, {3.818880515130136*^9, 3.818880556947299*^9}, {
   3.8188805880547743`*^9, 3.8188805966708937`*^9}, {3.818880652990777*^9, 
   3.818880698640092*^9}, {3.8188809110545692`*^9, 3.8188810293916597`*^9}, {
   3.8188810600490017`*^9, 3.818881077562565*^9}, {3.818881274486465*^9, 
   3.8188813326116257`*^9}, {3.818881435708609*^9, 3.818881460693994*^9}, {
   3.818881695629972*^9, 3.8188818143038263`*^9}, {3.818882259098098*^9, 
   3.8188822636145687`*^9}, {3.818882351819837*^9, 3.818882359491149*^9}, {
   3.818882392746964*^9, 3.818882397313724*^9}, 3.8212943095227833`*^9, 
   3.821294378858994*^9, 3.821307484977674*^9, 3.82130753209188*^9, {
   3.821307598578559*^9, 3.821307600938452*^9}, {3.8213076882876577`*^9, 
   3.821307694044342*^9}, 3.8213077810216227`*^9, {3.821307813015675*^9, 
   3.821307878802382*^9}, {3.821369625294866*^9, 3.821369702093603*^9}},
 CellLabel->
  "In[125]:=",ExpressionUUID->"3f07b098-bae9-47b3-8b14-ad5b359acfb3"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"dificitAngle", "[", 
  RowBox[{"11", ",", "31"}], "]"}]], "Input",
 CellChangeTimes->{{3.8211347464481916`*^9, 3.821134752542008*^9}, {
  3.821307498701044*^9, 3.8213074997655478`*^9}},
 CellLabel->
  "In[126]:=",ExpressionUUID->"5b85f204-51ae-41d9-bab1-f01ceae1173e"],

Cell[BoxData[
 RowBox[{"-", "0.834975480124432`"}]], "Output",
 CellChangeTimes->{
  3.8211347608356133`*^9, 3.8212943151462107`*^9, 3.821294383316283*^9, 
   3.821294616767941*^9, 3.821294828020338*^9, 3.821307500913349*^9, 
   3.821307536249859*^9, 3.821307606984208*^9, 3.8213076982557*^9, 
   3.821307728800078*^9, 3.82130778438398*^9, {3.821307842486602*^9, 
   3.82130785308265*^9}, 3.821307884762382*^9, 3.821308425008676*^9, 
   3.8213084988154163`*^9, {3.8213696580832148`*^9, 3.821369675829657*^9}, 
   3.8213697059226847`*^9},
 CellLabel->
  "Out[126]=",ExpressionUUID->"1775c23d-fba5-4c50-b85f-69933b890bf3"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"CreateFile", "[", 
   RowBox[{"FileNameJoin", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"NotebookDirectory", "[", "]"}], ",", "\"\<dificit12.wl\>\""}], 
     "}"}], "]"}], "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.818937698735478*^9, 3.818937756495035*^9}, 
   3.8194018484369392`*^9, {3.81965853232959*^9, 3.819658533362987*^9}, {
   3.819990708238576*^9, 3.819990708766407*^9}, {3.820323729422749*^9, 
   3.820323731590892*^9}, {3.820358665383013*^9, 3.82035866605772*^9}, {
   3.820404808908366*^9, 3.820404809491961*^9}, {3.820491107496114*^9, 
   3.820491108654986*^9}, {3.820508649384589*^9, 3.820508649871284*^9}, {
   3.820525520844355*^9, 3.8205255222712183`*^9}, {3.8206020615882893`*^9, 
   3.820602062069152*^9}, {3.820617865583745*^9, 3.820617865896768*^9}, {
   3.821308438670103*^9, 3.821308439067154*^9}},
 CellLabel->
  "In[114]:=",ExpressionUUID->"3188ce0c-ef67-443c-986a-8db52dd3481a"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"tempfile", "=", 
    RowBox[{"OpenWrite", "[", 
     RowBox[{"FileNameJoin", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"NotebookDirectory", "[", "]"}], ",", 
        "\"\<dificit12.wl\>\""}], "}"}], "]"}], "]"}]}], ";"}], 
  "   "}]], "Input",
 CellChangeTimes->{{3.8189377704882917`*^9, 3.818937854033463*^9}, {
   3.8189380304670687`*^9, 3.818938043118898*^9}, 3.8194018508889933`*^9, {
   3.819658535534758*^9, 3.819658535963414*^9}, {3.819990711456543*^9, 
   3.819990711826008*^9}, {3.8203237338332853`*^9, 3.820323734042939*^9}, {
   3.8203586678363047`*^9, 3.820358668093233*^9}, {3.820404811382395*^9, 
   3.8204048119559402`*^9}, {3.820491110635051*^9, 3.820491110902927*^9}, {
   3.82050865169508*^9, 3.820508652178035*^9}, {3.820525524264052*^9, 
   3.820525524464705*^9}, {3.8206020659296618`*^9, 3.820602066462071*^9}, {
   3.820617869277397*^9, 3.820617869564329*^9}, {3.8213084416078157`*^9, 
   3.821308441817216*^9}},
 CellLabel->
  "In[115]:=",ExpressionUUID->"bdd7f63d-1c86-4ce1-a293-c5201e0c6a74"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Clear", "[", 
   RowBox[{"show", ",", "pies"}], "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.819401896828451*^9, 3.819401908064958*^9}},
 CellLabel->
  "In[116]:=",ExpressionUUID->"d2252a81-aac4-4e59-9d1c-3f7bec4d2273"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"pies", "=", 
  RowBox[{"{", "show", "}"}]}], "\[IndentingNewLine]", 
 RowBox[{"Dynamic", "[", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"pies", ",", 
    RowBox[{"Joined", "\[Rule]", "True"}], ",", 
    RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}], "]"}]}], "Input",
 CellChangeTimes->{{3.81894215244585*^9, 3.818942189403638*^9}, {
  3.818942260424424*^9, 3.818942272990834*^9}, {3.8189424469328547`*^9, 
  3.818942448636724*^9}, {3.819658489532757*^9, 3.819658489757413*^9}, {
  3.819658587376926*^9, 3.8196585945181313`*^9}, {3.820323841752162*^9, 
  3.8203238509529333`*^9}},
 CellLabel->
  "In[117]:=",ExpressionUUID->"052d9d16-3cfe-4c04-b83b-ed3d9d4f54c0"],

Cell[BoxData[
 RowBox[{"{", "show", "}"}]], "Output",
 CellChangeTimes->{
  3.819658615528438*^9, 3.819658665260564*^9, 3.8199907591075974`*^9, {
   3.820323809688353*^9, 3.820323831802836*^9}, 3.820323863537983*^9, 
   3.820323934348301*^9, 3.8203446138239603`*^9, 3.820358728503972*^9, 
   3.8204048337991734`*^9, 3.820491126421528*^9, 3.820508664409593*^9, 
   3.820525536582492*^9, 3.820602113855201*^9, 3.820617921610487*^9, 
   3.821308499119285*^9},
 CellLabel->
  "Out[117]=",ExpressionUUID->"91f4d2c4-7ed9-4e5c-afe9-ceec0a3e1658"],

Cell[BoxData[
 DynamicBox[ToBoxes[
   ListPlot[$CellContext`pies, Joined -> True, PlotRange -> All], 
   StandardForm],
  ImageSizeCache->{378., {117., 124.}}]], "Output",
 CellChangeTimes->{
  3.819658615528438*^9, 3.819658665260564*^9, 3.8199907591075974`*^9, {
   3.820323809688353*^9, 3.820323831802836*^9}, 3.820323863537983*^9, 
   3.820323934348301*^9, 3.8203446138239603`*^9, 3.820358728503972*^9, 
   3.8204048337991734`*^9, 3.820491126421528*^9, 3.820508664409593*^9, 
   3.820525536582492*^9, 3.820602113855201*^9, 3.820617921610487*^9, 
   3.821308499130713*^9},
 CellLabel->"Out[118]=",ImageCache->GraphicsData["CompressedBitmap", "\<\
eJztfXnQXcV157Xee98nCQljFAfjOA7JOLNmaiZTM1VTNUn+SCqZmkoyqWT+
nBlNGGPMJsASZvMCCCRACIRYDEEsAoMECGSBVgRCaAftC0hCC9r3Xfq0L53+
9T2/vueu797ve5+wXe6qJ32v3729nHN+p0+f7j79t1ffdt33Blx92/Xfvfqq
/z7w6luvu/67g676q1sG2qzaV4LgK3cHQXDkqgB/G/un/OPSj/CPfKnw9w/d
fw0zc+ZM8+Mf/9gMGjTIPPPMM+bw4cNBD/db3SxYsMD85Cc/MXfffbd57733
esgrx48fN08++aS54447zCOPPGIOHDgQ8LetW7ea+++/39x5553mxRdfNGfP
ng3+wf1WM3PnzjXPPvus+1ZUVk1+u3DhgnnzzTfNXXfdZe655x6zevXqwvdy
mv13Lru3e3THjh2G6cyZM2bMmDEu//bbbzfjxo1Dc2sFPekEYfqG5L7N/fdN
s27dOjNw4ECzc+dO17/JkyebESNGuAbt3r3b3HbbbebQoUPmxIkT5qc//alZ
s2aN++35558377zzjvt79uzZZtiwYe5vlIHq8Bz+fuGFF8yECROKSJ5RlJC8
n/v+2GOPmXPnzplt27aZH/zgB+bkyZOFTchpdvA/c+g+ceJEM3z4cHP06FFz
6tQp89RTT4EMPaQJOT3qDIXixO9lNm/e7BjAdP78eXPDDTe4HqJVY8eOFSpd
Araa0aNHOym56aabzLFjx/w7t9xyi9m/f7/5/PPPzX333Sfv1MymTZucsOLb
30s5igNe6nLKk3K+6qizdOlS/84TTzxh5s+fX9iWePtrvv349rfu37YUJ5Yv
X2727t3bQ34GaR544AH/s3SPcqN6l6itVxlqBX1CZtwpxVk8meuuu86JAfBl
FZHH14oVKxwR9u3bB7TUJBuNQ6PmzZsHmrKrYN/3v/99VCZkb0uRPVZSmy+J
CSC3lGH90DrvvvtuXv3SXNbP1kakrmlS276GlUJSoceYYt1o6G6UI0hNd8OK
mknKOuRmz549Ds1vv/22Zy4QZ+XJl44iBg8ebLZs2QIOs7sUw/fffx/MZXcB
ymuvvRbaLpfcsZLaUgJ98803QyRYP0j9+uuv59UvzWX9bG1E7rSOmTNnjnn5
5ZcdTaGVodMsUymu0iP/uOpSOeLUdJeSpAfJIfRDhgwxt956q6sceUio9cMP
P4wxFr8DwlbZakkDHIH6hKQDM7a9uaSPldTmS2ICghOSPmnSpLz6pbla0tHa
ZqTn2IjhA0VYNU/SS4/846pL5YhT012Kk77mJB36BwSqu7yvOul69dVXfY1T
p051NUEBQfsfOXLE5eNdEUzHdDu4eG26ceNG1yF8+ztpHfQhOogyPI/ai4qN
KfhFixb5Jo0cOdLZFEVNkm74JrEb+PbX0qRRo0aZAQMGODXimVxziufxxx93
3xrCB+mhb4LqYaKqXmUoRu0empVXOjHDaI0BHh8MDEjAAYi2a9cumJquSipi
jPRQA6gGA4roVKcQoaAXL17syoGZagcgQUDN5YHBlgnedswoyxs5kNCHHnrI
iR7kSBs5eW3IabfX90AAKG//9xwC1/Cc7b/Nw7+/7UjGlNOvztCpr6Z+w4EG
+owfK18eVNAXKA24nTJlCkcf8BSqD5XCdrUNCPTID9DBcMXIhK7QwsRIjw4Y
ZWFmlUXGgABvvPGGM+OgJlatWlX4Xk6zQxC2mSVLljgLQtP0+uuvr6HbVnKn
TZtmdMrqSScII+R2RWbMrX5LaAPSSwp/v6Y1v/cN8+zPN+GfHvlv9fuS8i7/
JWpLV/LCMeQ3pP0NaX+V8gpI+x8lb+f+DvO9odPML2atC/5rxnM3ZOTdmJH3
zRJV7Ttm+t832bz14efBEck7ceqsOXn6rPv2J82L2LzriCvin36xPPhTyVu7
5YB5Z/Z6c+joyeBvJG/j9kPmix2H3berVFUDH//QvPH+muDfSt6Zs+fN4WOn
zOmz5zzTbxj2vrn/xfnu2x9J3oZth9y7s5Zto87PaiaLtcXZXp2Tb+G/i9fs
ck23n+BSyRv+2kLz4sSVsarQFTxnm++7c8dTH5m7fjbLfbtS8pau3W3Wbzuo
+Fwzx0+ece++8O5KLw9ox5pN+82+wyd8mxd+FrZl3srtmd3+luStXL/XPP76
YrNl9xHpds2VhXenffxFESm+Lnn//4Gp5tExzsL1VaFIFNGMPAOGf5DiRDPy
XCZ5M5dsdXWv2rDPWgfptoQkazPHTpx25Y0YG1rhU+ZvBJf7ys+WaubzLQfN
6TPnzBYRPtsElnjKZo+dvhoElRlIBnF74Ovvyq8zFm9xv366MWqXxR5AEWPl
8nV73XPjP1rX3WT+7pCp5kfPznbfCBZIJaQT33oJLW57YqZ7d/7KcJK3bc9R
s27rQTbrsy/2A1tmvxW0/+T7/bVEv61IFfJj76HjbKCrRGSA/ACZ8BN0ieIH
51zb9x5z/BDsm6PHQ+6+NGmV58+U+V+4vAmz1vuuzV+1w3z/ofcs1feYtZsP
uDL2HDxe890gTdFny1gUHPTMYCmtWwgGGnHIKheIpqWddLHufnvizSXunYEj
P9TcdcztI+WOmrDCvDrtM2limId3bnlsRozPIFFTPrte/L6UBVDhLeiG2cu3
meuHTTcfLd3aJj+DL3j5Q9sxpEdeXYg8Mnri3A3uZ6sUhYA1R2jgBHRpSN7P
p35mZixyRQTf9m2w5VuZdvKNBGKC2ucvXGDbD3eccmJhW+TLGj9znatzwaod
IS1cWXSo5QgNnRgb7Ggw6p0V5uNPd6IednPMe6u9LBFuD73yMWUJChZZIFCe
LD01bqnLe+S1hUKKdvc88hav2Q0l6f5+ZcqnEIVwiuGafpkSFAAVTCVRKVg/
fW6uyxs3Y60XLAgV8sCEpGCBASRDQrASaoMiduuIGa4WLToYJ6EFX7fjpGtq
X/UDBkLLnyBCBie5YOIc23UreeFrLI/6cMGnO8KeKVGsOcUDBTR2utRm8ybN
2+hbSr2QI5P95KEDR046AAsJe0r2otW7zMg3FrP+RiSFNacrWMsVihuAOHiN
hFHGgc7+9vLkT13eoJGhGsTIwj7uO3TC5VkieJFFz5D32rTVKZGtmx1iDInI
KorVHfX5m4Vo9LLwePayUMIs3Tg5hbBpYd19oCMkZ76k1h2dkWfHO0gqCbPs
8z2eMKrNNQczsNmaAFbOwzxYZMAzXmUBEBIWQMVFarz5wVpPNSXL4bsc/JQg
B+y4ZZ8rEQpj0twNPZsLNyCVFG7SyZfcR2gBGtz3/Dz3o1XoafnGD9DaVeSb
bSe4Bczhm+1uAEODYIVSvmcu2RIRntoVdh36b0VTCB8Hu32JEk+dkSXxFI5I
4uvanAPj6+7RS8zoSavcc2yWtSgd6wma8gJvf1hiNWkkUMYjMyn/1LzQbTQ2
BAPmXsub0ZNXhfqAdEXBGWDAwIOskW8scZVarHnZhBIDHTbtPOz1rMYGBwdM
SjAzOHf+ghkyeoH73Q7ItTyMhIU5xQFaffLZTvedXLKQJmI4JmjE0KTSiPnh
kx+556wJ4MV3qa3WWauvL44jBkNbHmLsmIYpXgwx7RoxLsFSQdNAGFa3auM+
9kQNXawzhSWZbUA2rL3s1QxUNVlZBVoNN2j3D2eM5hdWydhhNIhpQyToa8wK
7EzWGdkNEeHrHp4e66Adt11h4OrV909Jo42sx+wNNgGmBii8f2iGZwHvvY83
eXf/ufPurSMd4dABcENXYYAmdcTIjwEBAsnirGCS+5iXQ5CsVe39yGJwu2+X
yOvEL8YMy76G1DRURFbzhjIL9R3T6znAbTjVppPleEAEIEU1hDSgJcpPx4kz
zrlgJaAUcO203mUNffljV74G7u4DoY33hh1AioBrf3d54EUGcCfO2aBYEAFX
KTT3vSRwrQD44hTp2mF/xUgHQxgvW6kMYsIPkwMiaeU7qAt5UBNSOeRiuJWm
euSC/7DqTSZyJ8/bmOoJZLV/aLZmItcKY4TcUOebC5EMJMBbl1LhVSB9YhKA
pKoU2u09eJxjHhQvXC7AvMtrizRDew6Ytyq2Rvy6Z9TcgKugpBUokIFtcgmf
6Z9solwq2XMNhwerIXSAEQXwbd0dLj5YS9S34Jnxy+K9V/6T9xduDrjhgnTG
NMEaSOYfB09BS+tSPeykNXZWyjR3xfZUJ0NJjgN9dRzobRroTrvBar6AWafi
q1AxVlkIHscjTx5aCBbzxBa8bvjJMiTObIAOroEWwj2kTU/MG5z/B13qL/Oi
B16an4Q9WxvJR8M5DjXsLdp89Xc/M9s/TzWZpVPz8W9FF2KBQRIphX9b7C6r
45GsKLjfhv38E1EBDTdr0eTuOHlGTL9SeoDaBfbIj/9pjqs2YjO1BJ1fYV/S
ekD7FWFYwmXiSq4JY1E0k62WHYM1Uk039BB23vn0LPPM28v8m0hJNZHw+FE6
oDDw2NNvLU0qDE1wrS2sNes5KZMYN16TPqSoxbh/7idRlwNO5iDfFdUGh8qP
loZw02qDeDBq/ku1YZWZcwqdv3ABaoOvUB/CfUDLAtxFD0AoS2tm05Uine0t
JQBISKHB4ZVGr5jSsLR1f39qDUMQDOpXQ+pnb3mF55Jlp58/a2n+weMfSmv6
e/uXguh63BCFBtcAU4a6YPc5WttR1++ygjzABIaTXGkOTrJhd8EPhuwhLy1g
ScyyuiMoUh6gBBQIe2SJ7Pk34NEPEiJScz1uqjwI/kfHLHSTQioGrTRAeiT+
9nBCabw983PfKDhbmygNWmF27qugnVYV9H5pVUFJwlJAPfEcliDiquLZSFVY
0dF0zFcVTgz8uvyK7c4kxmPOT+elxYuGVeOhiujpVQSAAlrif3j5nCvu6EmK
AbUGPo+OWURXGsxvO8ni1DWtSqguqErs4OD5yvlnfxrKarL05LglIf1ocmTq
jrr3jcB8tu1qy1EbmK2xSdDQEHjQCS7FFZb7VEW6Fnzs9NuzEX5+5MGhRYuE
Tg18qKetiuIrWpppdRAlXoGEjmcUC4eNMzYwSaI+ASPoBqdbH5MmE9kN4ICs
o8G3SDHXmgQezNZpkktBtNjgs+fgcW/eYTq8Un6HaZfUKm2VdEjNyf7qUK8H
nLjBFORzxFRad7Rr3RFTBy9OlG00bMxr4uXO0iFgAZWX1iGcQ2XqkEZKh0xb
8IXvk7VhyB+FYzJ03Iy1jmoY4q2WbJeXnpuwwomHHbG85YJ5NOq3j6cVCcXc
9pZ9gsVGG1MGyjxFwsFSRoywhN6u2cNe/YSsc/5kSAoWh+jCmrVsGydpQDS1
gJ1KhEDv47wHGmkQebgZaMDTppi64ItszfLIa9GkCzCHZok1Hy4Aq0bIVUrY
/sMnnDkbSU5cuVgJ8ds3IyG3A3z4uh1H3HBuS/cT4+R8XlYJnS612GTNymAL
yCUtn1aZaK/Bz8LHvcgUKRMsrHLZSL3rPhj0bYvIjzdnrNU6hDVSe4UPpRWI
AN0rEGIiU4HYgcM9cP3D02MKxEp9lgKBTEKBk5WgJ8akZgqEjug8BcI8UR6q
bzUntc0VSK+UAiEoAC/IAeanXodMS+qQcCEe6fanPgo4OkF3aA7Z8ayK7gB5
IFpZuqPmnCHoHPDQW7oMDcqkzEEkWd3zhUN20TfY972z9EJN6wUhca+UUsD/
8CvDpQ82hb0KJSShEdy7nANatUBDmB4/aD3RCIFueYyjUSvtr5Ch2K9UCTB2
SHRLI8oH+jvQm9qlVcIlUgG9zfioaT7nxlgHpNT98MmZCdLFVYCdOcZUwFvl
VQC8CVoF0GbIUAE3PvK+yxafdVUV4IYjhflwa4S4iVKzM0tY72rXoLf056PE
vPYBU5gw2c8EfUTC5qAP29mun3O2FuThoNVVliDS4565gAc7rQhyLKZtpLH+
pkxgk1jXBITzWGM93HlbI8gVCet+54Z8MmGOughzeAVyYC5U89S52eJS1o78
0AYzklSMAN+WgfDnBOGgSRbCsdUDblU8AycHLT2Y2wrh3IclguX0YCcRLm/I
VMsDgWvd/AjSGyLHGuRgFFQOEqYEfaVYhW5vfaJ/sK0i6tWdtsOKGQ5m2AFL
w/jaB9/rNIx167EYaA15ImGOrFNqvmqyZXkUqmH4QcHwPw6e4vO0BmStStd4
5zzW5TApwUzAThO7At82L4y2XLgg8QNsOCtjrC4BWxZDM1Jzj0KhoXvzox84
TOAbXT1xPJjk8ExywtMV+s7zUXtZ1INGck+D8ylcptoWw6DwVOYtkldPopEY
wOayOHXrfjjGxyq9dilHL1qI8ukC+ACswS/Mj0mPAl/YmCzY4QNvRRJ2GJti
/E+IGZkEx4TYh5nNL4M5rTGSYMEwYc3YrmLuAjFnKUaDRpIGoCz3x/zrmBnp
PI48ugn3vRChiuURKFgaqIo+LUcx9NGxBJ+tRZ+yCxtuGVP3GaDgcKkx52Q9
Zu/WNKC8jAOTBUBKbQ76WgGGOuxQg81yeRiyaKAoY1WpCEOWnr0Vl8fNiJwA
azcf8KR+bKyrza3WZCGHHm6wWFPNtl5+a3OTBzgh4GiAvsPczQ4A2m5FgpeK
uOG2OnzrVxIjfLcsHux41lU82CZ5GKJ3cMrZkTbQYFiwakesXA1/ypUW6HZV
HvPqGc9VF3wSDc4J2erjS8YCOvgKLU5Cas13aRNp71VKti+XkrXPlLTCKB3v
b1y2bYtpIWo5tmpIy7Ga0LudZuSOlllSadaybfH2Z2l7Wx0mk5zpWTmm0GAt
/xeyGa5PgdxC8yXlVmu/fLlNCxy6RAHprNAuswJKsogyQbszhfYKyVPw7gZh
lFU6XzKYJYwMvl1N8AIsJnc4QyNcc0/T5fekLpi6pIkXrpr2yMncOCShktmY
ZD0XCWOhZIn5HOuokiqvcbH2N76iWPXLYbVmAzdya8x1VoQmqf4TnzOXbPHl
QrtQhL4teVjLTNIJa2eYrYSbjHPEReoXA8J9MCb/q8qCwfmb/b9AMOJD55Wq
BSNejybBdzw9K9SJ4E5Jbg+ILIEYt1+cuKpT3FardKW5LdStxG07fnvLPouz
aheq7yt8NVgL6q+UXpqzYf1YsMO3f1OGo0G4URFp/baDweVS67UPTtN0+P0m
rByu3O4wVbJYeaXkYUcCFjo0i+R7UxZxCZ0f+0whi7DYnwSJKiOXRYOemOlZ
pPV3FotostEI5gdLeQQ1kqxtlGYLdsxQmXw9Q7z+m8oT+vs2a5En3bFMCP8Q
vv1BBRpboz6i8Zw0jUeISRu1PU3Pm4a/H6Mnnyc9sRUmSU+qlzx6doR+ootC
Ty1bfyZ50z/Z5POukDZz0wo/ll6ermquWY4mz85J9SEP8mVpoefHzPstyRuo
5sJZtPh7yYNPJdlvLrjrfv+PjH5cXqbNrslsKVZT+PY/SB7W10U6ghulBSs3
7K3UAiBY9n2XpprathB8VKZeJdXUamrQC/51k3qXSB1cnm7Wt67UwYN9/GBA
Zh3aS/JnpcqruRkh84aXfEebft3xjsP6/+2+h+QkTeInniXGRh0pI3y/f8ZD
mM+35iEcNmn6EPypuQ+Jex8xsBjCj/FrkKZPn47gTH6XTFYMRG7ggd7Kr+cy
F3AE0QMRnAV/IwQPg7ggTB2i6CHlBP7TRy7yqwljh6ErTAjXY5scru61F0W2
8zsriqtod0Hz0HomBMNisCnuhVdB23y50CMSaMbmUY7+T7xwBGJ55ZVXYg3E
o2ggN2nnhCwjibCLq4hEiN/y3HPPxdpvacCdlDlRucqWvnr1anPvvff60tev
Xx8xIB54qq1kkaABgs0gjBji/zz99NMQi4ZIZUYIJS+VWLQoKpoiApGB1L30
0kvm9OnTLDorPlCpovv90j5U9w/xk6fNVpapkUL9a/pQGXquLEvPMjViBvjr
+1ApkYLt15TofKgZ0TGta1ojvGe/kg/Ba9aKh8LdemXoibMzrdFVWKH9lXqI
wy8ismG8YRIrqSiyc618+X1T5WvrKCfob5UKaroCf35CYtj6SiWIbbLcfNsl
v/B2BnD1hetAtD0SNVRseXtRjNKuF14QhTMqfHtB4eIpQRkMOGnr8IcwcqJL
1soU7Y+9Wtm48cYbY+XmhEwsVS4fQuDNjo4OEevw36y4gKUKbTVW61okCzXW
xVYhpajQuocuZu9+Q/Qv4aG6f6gp0VvWrI0te+iXkZ5lJ0C0iQprLOWpwVHD
2EOuOvsDdoNL4lYzvoNdNS175/l3VlR+B5tOEu+UNIGK/ETVqyu2iHL8RV2o
L9sSKOM56mpdPKqb5T1i2XEL7H93pjNNPUi5/RhXrZomrqQq1RTaUdlupSrF
Z9lSbfJAkX8JKxxVKZVlVxU5mqrU0a/ZOxmC0tJ3/kLe0ceErmhWTt2Xo1Xy
n/M1ta2raVEFTfoVLwe7PlpRDvYOsJxvyeM89O/L/k7ryh6qdql1pezzMq/T
ZTPaaWfLzlr4+Yb7t+7j2OHzLqLi2D/+ZeuK12dIbfHfkFJ4RMf36JJOVOlq
ZEUTZq3Pqkivba+UPGwoxOlJfPsvGc9V7Cq3FOC4PvO4+4SxUBKt0nE5x0qe
OrwVrnXmtBQbZrvYUoYqxbdwV3HDyYE+gY+l2azm/qHkMSQcvt0heeqgqN9U
oWVXdwvxGL6Ebr2ruqWY47t1zZBpqW6pXaW+W++q/SO6W9hjiG9/JHl6gOA+
B70TgHlDVFQZf+hArbpzb2c3k2RntC3Xk+Sm4RFW/5fk6S2e3FKl9yVXJZPe
kECSmKiNhWTSW3kqkokb0HC8X0gl0C0GhCbT70meiu7oya9J54/2q22kfy15
WnWQnCvVQSS+q8M8aXIiPjW+cX+H3oXDnTl3PD3L5/XJKK8q2bE/N0X2DyqR
HZu2qpBdTgXHyC57JJuSXW+g+13Ju2ZoBPZ/J3kq6mpwleRhNQBnvTRJdHBS
zQpuBirLCrZZ69MiVmTlDc1gjwoo2CX2hAGrfFU44shNyNwjV8Qy3VyybJk6
cljEMkZ8wrevSp46NBPQZpij9gh+XbUzKe3a/GO92AcqbQyR56ZIJJiMvb8U
TFa7MruVyY4APDkB1hrN2rnZrOVWP2zzRwgP3UhETC7DWoTlKsNavW2PzyH8
nhhQvl7s/M1n7V6r8suwVhsEmrVLwrjgnrW6R61mLXmlN+smD2Un2e24Xcxk
Nl0ONqrqqzMe6pDEZhe0KvVBQsYvc/hHoDMeaNOM7yklKKSE8VCF8spIvjSj
OZQHbEKYKsfN2Rycr8qXByRRYgXyUNfygBIpvHfI2XTdVz2CtwldNU3xCePr
5YuGPhBVJBo8bFFGNHI0AbbS5gsJtifeK6fsqwrJt+R5uGdh6mBzMRLOWxQJ
C4PRoc0UFoTkknmGFxa9dVUJC0PBJIWFO26LhIWRtVonLLaD3S0senM3n8OV
JUlOq0MO3SMssBCSwsJQuNyD20xYMIzohAiw7BQ84DKL8gTklQDGeZjDPH1o
j0TKEiAZspICpLYriwCFJw95OG7mki0cJjbtPFwoQxw6ECQuVmQr5Yp4uRhy
JWdzYkKhja6mctWu5cp9MIGEiW/lKrRmuNDCwYA3QljR8ocrkmI1qQtihSRn
pzz9JJZWTKy0TyNLrJIh9SJ61GO0V6KmQg7HRG305FXumKkdIylq0J4UNdaE
xTJAUmtKJDG1PTcYsQrfehWKH2VKR0OiIa0dInxXiylr0+YuT6Vr0bXTblYD
qHW36OpAQDxaxrxKoltPii7XVTQXY7EDGDqLodSrSLA0zdvO2lXpDkpwdShL
1SPhvBaCc2aJdHsiT7M0S6Rh9spUuplIa3Y4keZKT9reyxbqUWo1Vwu1nNHx
fdQzrWKhpqSoayV8C7TnmYKOcI9Zgi4WWGb0jSxBt9LGg0/alcKOUvgt0LtX
+B1h+AaOy0mAz+4QeUQ4gbmO8Lh5Is+j/XZULBR5cRa7b1SD2tnsOtwmDdVT
B0odqQrYJ0UfYZrKij7uaNMRv0qIPvtwnVK9bREcGq4uzVVEI2KcA4x4OP8C
3W7hQdZdXGQwLKtGBiOmR0Rrd9xUZoHnMsKPtWWgwuolokK7yUlbBIuBoSfI
qAIMy6ouAIPhq7sGDH1iuBHJZ42GsZKtCBB0yNrptF8FgTwjIfw4O6ABQcJo
Yze380UgIeFV4GotLlVwQlJl4mRZ63GyZfcR394SOMFIqUN1F+DE/UWdDgOb
1S9es4tUQLBtUqEzAIFG6QRAFKG7EyA6JB+J2VqoqEO3XYMKQyMhlia7ogsn
iRAuuggqlEV6h8PWN0cOPAkZyGHw6m5GjnYPKuS0pZCDiamVRR/RTYFnWgQe
eGRi4EkGJcfHIucSzUq9ChbFJiV+1HKkxo/cvxTHz+ETiQ7W/dIZ9lrYEZEk
okQgxOGz45fzaRU0TCNJnxhvgqQVXUASA/Uo1eI+cHtQaylE6ciqBBUuirjI
oEJwP1ijuEX5txXVmSywKIHKZ5KJK+3xJXP5bhdxRQ5jdsuLtmyivZIFLX0Z
YSegpWe8hdBCmjhnQxNozVFX8H0t0du80agiuHYzW92HiacZVCoLX8MVvhT1
E5dYNcFXGgRfGr6USzALX5g+dR5fKQLm46umSRioBbb40g9DrSYAx57CbQ3X
BduCJRTKpzbbi4An9yxLz0yngOfeo+cU8s+1W2s3kpzYIBH1t1NwGyThe1sM
N1JH7wfgcxJUJsrLhSW7zzVGJNvnyBOekAPSQGPUSjZ942ABF0AEp4EGKn3p
vPaBjmnYQIwbIaANmqBWBSgCaqkzeUlmJECl7EsyDOHt+TpDK+J13sGFj1xY
zT6jT3Sk4zOJwA153e4gqaNHVoFtGHa67mVEUrvGqTUV0jh9PAeneg6chVM7
IcvDaceJM55OWZsCeBehlpPX1RUWCz7dEZ6ziEPVjqGdh+p9z8/7kqGq/fAK
qtiJxdZEy4CrU1DVmwqqQ9VI3G9CVcdeIQ3UFTIaqpBI3A2PxQrctWE8VPs4
qKoAju5vXh+KoJAvTlypocqa8FiSc0H8fhKFVN4eWgapG8ojVQk4b79klwk8
CSAdIZUHIFQUav+wpU31sdUxMQZShufEvZjsFu/FNEmQascyK9c7FWkXZAE3
FnsSiXEmWauWKT6nMOd5oa/tzoHu0NELWgJdWfERQnQKujLRT0BXBWfNgK5c
9+aSgu6+w9GyAommrxjlc/qS4mLospTvqX161LL6zg+kxP04oe9Hovx54qj9
XBGmQ29HTKD6SiVJQFsBC2lpy8oBsyoLiV0Y8lLo9reg9tKj7gjyBloX0Y3b
Hvh6a9BdT6I7BpTOQ7vdQ5t3otGhr++R06HctNs6CUQVpboaymV/TWWU60hm
VqC7DeWqhHiwagnkCEujyIMFFKwKt9b71ZdrhsZGQsBd+hVWoePXY8/k4Y5T
BJVc4REjT1ncb44uDvBkVLHMNWeylIAyqbwS0BeeIHklUI+wF0e9FbbKqN+Z
gfoH46if1DXUY3+ecM+ThsFMkRKo59AU6TuyD1iHDkA4NGw1wbYUE2GdE+Jc
rOPdJNZtpbfLtjsL9/hNTiVgDosOl4lwvoxbXHA5A654JMyzQp9raSbt9DGB
LsMc8+NmMEfiMQzLkCTMEWvQXfhiKsAcO0xKw5y3bdieN4O56GkPc8xDdIpg
HnpZNMzPnD0fcHk+C+Noc5w28ZXfIoyjryoVYVyfg+BcN4lxaElQZdDImeF7
Mazvibb1aqxjvYzt7yrWOcZKKsK63iQVTTgouTvl4mekBMKFxymEx9cIalFp
Db+1JIXthpYtfTWvQnVZUPPoNuxOMGHfoRP+IHJZEKvQ/xcPxDSAtu89xuCF
BLFVikHsnvKuIthV4NQGxiyFYLfhNSvANiELGOsjg+P9rKfhNo7EcFASzvoM
A6nEpVyj4EyDIg/OfBdHN/lujEVZgEaSzTExJR6T4CBc61f3PMQgrbY7A9Is
WlkxYTFdgDQvxpJUBGm9fJSGNM61MXUe0hzNcayEPypYa/mPvdcikH9NSsRe
TOmut1ReSt1FETmVI0nuAqAxG0Bc+bKAttqEgGacbdwXUw7QvCwKqgsyoECt
rp4Meulb4Nz/+I5xAF5k+EP1zQGaRFzV4EYE0tuY6P7foRGwJUBDr7wZuPsb
UnHBYyd+Bxh7rACOzrGqshhXQ7aPGcEdkrhiQV/dlYQ3FuX4Lu27uMoPf1Hu
cdCeJj1ng/xY3oUKTsNbbcDtNLzbPLxhYcDFAPe/lfeeCYA23CWiHYqisX16
hdBGuMhIek00Qqj1SJf0qM2X9YaLXHxzG4RtR1l894tEv11Lq/tgiiC+Wb8n
RYN5256j5cEcyVAaiFbGtM8KEi/nLIlncvmLHYfjeE4QL8JzmxAIJ5pUPcS0
gnQjBmnIzcg3Fgf8DWMmLpYW5e9yqStGTSiA8ssRx0Lc9CyE8oBHP6hLldjl
NyhaxcOalRdrfdephTB7uTl+E5+7v8BOUgg6Li1GGy/jmyurAlbfdV8GsKxT
n06loQMQ6xezQDy0EMRhdk+3Zp0FZOwctDOZaA8Xway9l0jlwaxXSiqDWTn5
88FMr9KnnGoinM748eN9SCRep3rixIlU2PKsiJzkXZYu4Giq1o4T2qGht+W7
z5rNBwJ1bWpcNdgBMVIN9jeNMKsdOOCUVQyXCK1y1QCPkUmK1ECQtc26jFao
uTAi2GBvnPoL+3Hy9FlnbogpllQEXO8zTRTBjmJF0JDaDkXRIaAEiIrYsgJ1
QCOlAyByPA8A5yyYe98LVAChb0GF83AM3nvoOOdwzXBvSdGNuF/VBdxjZza2
3CJd0Hs3m+O+4cWPNg7sv1qihK7gvy16Og54sAJ02rTzcAT4JUuWuFhiYUhd
/NvbPVwiQm4W3q99cFoM70hJvNO8B28iuQ3zEOiBeVmwP9Qa2H+waHMM9vCO
oDnW8NKwbyRh73bPKidmzXkXeY1hGexfhsxuRnWHugvaojqmx5RVDglkbaCx
wNu3QPmZPcz10V/NYwVzum2y4I1FUy6wnTl7Pg/eW3YdKQdvbMjrKrzrSXhT
PaswE6kBuisQJ32hZbFib2gLeYhvLw/xmoZ4iF2L67qL3jd06FAzZswY/B3q
3NyY1FlQttPNclCOLalmAZhnU5HI9egK6jVhKDBTCcAr1u912ZdHdTfctged
cB9ZDKyYW4nnp/Ngxegsf7cCrBLrpQpY4QjoRrDydZzzFKAGeUjFNWd8Pr1y
Ho9UxyZw55v4+TnIag81Ugsw3PAYxk2I2BkD/Z0cbNNYxjXPSSzDeIEbYQdd
oRGW26OG5IN3Ymnwhi5WXCGD0JwIZb5hAw2t4sDv4qO3KG4RchFoJElssltH
i9BuWZbRFLZBeIZCXFJhO+ueQvVq+Kx5fJqkr7gZOOkTqwTOnjFwQrSwhwt7
QF6zSsaik4ZcETBvHTGjM8DEljDcVXjeRcwPf1kaxTlKARNA412MRcDkJF4c
ei0E5mgPTMb+xVcGSUGS87gwg4I4Njvi2JyxeEtqCu05WwWY3OCOCFIRMJVr
12Iz6OPyvPtz+fLlLirv0qVLgzKXJdAewXfIEpTI2i0HmK19MGx9v6jzNfPk
uCWKS2Eb+JKVsfj0ePk671vx/S4C610/qwZWLKajQcb7buo6wFwKudiPy3TR
kXshCt+WQG7SM20cMLPmtvQNdhKuPaM6G445+AHnhiAJUVTldN22235tA2M9
pgaCrouKSZ7qwdDVeUyS6jwi4MWtFDrfTaKT1zQqdFpCIoi0BIB2ae3atbj0
o8ydI8SVi5Es7aN4yy2v+QDVB29ZEOMvYL1aHZoJT5umSVMEUbrEPM1sVzVE
GTPAtvNKabK+ez6JSD2T7h5EIl4akkJkzN7RpyNiBWdhUoVNS2Cy7jGJle+n
31qqvfvlYVnzsLRFSKj4AjjGhv+uIRBqIAOBcDccUsZoDgJ7pRCIpVNsAOGM
0xokQUznKUdGCfTxbB0WR2PoY9zV8TPXefS1m4ULF5phw4ZFrbGTT9ik0RCY
ey0Pp2Zp9DWcJaXQBz9L2HaGOFcxNYHEmHOVL18EFPL+YkscohBqW+3P7HYU
RoF+h/38E5+nUBiUg6EOtEo2qBDCuTCUeUwTGO6XbXkYSGMwxLI2xBc3k+fB
cPTkVSkYXtIMhtyzE+1SYRPhlFAI9PdEnDvvjtNb3Ab0Uaj72wE8igqPisDU
jtpE+cNuYlwzbIzRtK+ANWxkxA8LP9ulRrqDBw+66xpwKSguzMDtpvfccw+P
R+TcUlUEM70IgFQJZpD9NZv2i5MgzGPcVNt9j6EU6oIokmnUmDaNP+dQuGfU
XCwaKiFNYw/8eXPG2t/xz8TgBlXL1Dq4SeguzJm6ADfMG/FGGJzQtBhuuo1x
uCkbLQ9udoJZAm48GsXtA2Xh1lc6pKHVceKMX02asdZxf8/BUtAKouPKGPj4
HJE1YVYpZCGtXLnSDB482F2t8vDDD5sdO3YU3dNWBVZYGVOwavMhFDBmS4RX
eJHJPc7GB78wvxBFmARx8YDtIYrC1bgczn/b/xATViVjvo7UXiUHHI6ocK7h
VwuguNGRBRx4vNFcO0K2ADgIW0CCwwaCOvh/FxM46Z0R8MwMGP5BzPHMvcYW
2wGPnDIy9WzpI7YlYvnTdiMOHYzJGdDBCiFdqXZG2UeqwmoYTrdomsMHJZa5
BwzPJONkNwHD0GL+ABwfhk7CpQl29CNgdMDOpOhzFyWojmSN+7jo8/A/GqH2
rCZFX+90uxiiz204n3y2s6To89dC0acfobqYwxoBn61Aesaxyki8e/pNOyAc
EkxPrEyKfPfQwg2YgObWGKsm3BJKOSXc5BYaSUlSwi2NbNNjg5bnS+UVnou1
shxkCTPXarMElwcvpy6QyINfbyak7d63S7/EZ1/s9/qWcjpbQoKjOxlyCkMR
UzQrzoVyys2trZJTxv/G+FUsp9wrDEf7ADHLrKy2WlR7FcimamyubPYrlsPk
xlHjtIw62+k+1h731MWyIkwCnIenkIBTAMrqTfvLyiQmLOA7PO2USaw5HwqP
ZXiZxHCEpbHTZ8/56rhZ0lrHXjjUMRmRsUIRxYImNtVg3RJJj3pZIorr2UDu
rbuP9paC5qtVCiuabDA3f1n7rbL04bAYNK/fn9VU+pBI3TLSR5GxlpFXe9wZ
hL0BUQFVBA+bK3HeqNOCpxbVajzZ6h67XHhBnF9Q/bXS9MdSFoOplpEIJxBs
RJrvYRB8qEzQgyYB7Vp8RoxdRGYziIa1wasxm6YvplAogFEC7n5mdpBiMMeg
16atbspgXmgyZf7GTjM4vJkUTdq087AHPk+gQBCrM9cWwzMMwJ2VkivkUXqC
x81YG/yxsBtbsbAbYO+h437+MWvZtoKxhxUD2KDBPjsG/wupGOMPXUYT525o
GfO4mY9HoLG4nWIetYg7xtYJ5gFrUOrxCVsW86hauVUZzekEo6IoKZjxfUPy
lDIM/rPkwUM8PHS5BQWs+Q8UKjnhAJcreQB7FqMHLJzW8sD28IdP+msm3Kcq
O9A4eGBgD/WVPLwrFyZ4Gn4pTPib5gTHjlcQ/eDRk165jRZFg3NnVxXQFeMw
3sVcG0DBHjDSVUXZ9nTFEWU0efveY7l0JeIxseDehCIa4ow1hAP07xQNYWAA
P8dPnvE0RHAp7EXefaCjiIZ/KXl0y8MVmUWvb7rn6no9Ge6HmPmfRQMEoOGS
L7vAMEtQSN9p1l/0g/0d0LwfOOOMOM7Qi4l+BNGdSvj2TcmjBGNozO3OH/ji
QxgijOd3JI+zF+zr6S15dMFZMSxq9XWpFj42dhHjDTZv4eR5G0u1UDM32ULY
b5VaOCKjhZgJAT2IOBJr4Wxl7LOFkFZEFrOQ8y2EfwFOOay1sIXYRIK4KXZA
rNZCLG1JhPjgD3VrYMrytflNSgzCuI9yWjD495KHcfd0uEsw+FgV4UrIevFU
83rQecQfb/K4+/eaX6Pfo58l1/39o/Dv4Cv/DExPJ+c=\
\>"],ExpressionUUID->"4f6ad3a5-47c7-4c64-835e-8b8f959bf087"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Monitor", "[", 
   RowBox[{
    RowBox[{"Do", "[", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"{", 
         RowBox[{"random3", ",", "random4"}], "}"}], "=", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"RandomVariate", "[", 
           RowBox[{"UniformDistribution", "[", 
            RowBox[{"{", 
             RowBox[{"5", ",", "20"}], "}"}], "]"}], "]"}], ",", 
          RowBox[{"RandomVariate", "[", 
           RowBox[{"UniformDistribution", "[", 
            RowBox[{"{", 
             RowBox[{"20", ",", "40"}], "}"}], "]"}], "]"}]}], "}"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"angle", "=", 
        RowBox[{"dificitAngle", "[", 
         RowBox[{"random3", ",", "random4"}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"show", "=", 
        RowBox[{"Abs", "[", "angle", "]"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"AppendTo", "[", 
        RowBox[{"pies", ",", "show"}], "]"}], ";", "\[IndentingNewLine]", 
       RowBox[{"Write", "[", 
        RowBox[{"tempfile", ",", 
         RowBox[{"{", 
          RowBox[{"random3", ",", "random4", ",", 
           RowBox[{"Abs", "[", "angle", "]"}]}], "}"}]}], "]"}], ";"}], 
      "\[IndentingNewLine]", ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "500000"}], "}"}]}], "]"}], ",", "i"}], 
   "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.818875191515671*^9, 3.818875204914426*^9}, {
   3.8188752682532253`*^9, 3.8188752732118073`*^9}, {3.8188756772599993`*^9, 
   3.818875677674365*^9}, {3.818875861543709*^9, 3.818875869047735*^9}, {
   3.818881892265394*^9, 3.81888197347786*^9}, {3.818882003877234*^9, 
   3.8188820081295633`*^9}, {3.818882139855406*^9, 3.818882160002954*^9}, {
   3.818882298145385*^9, 3.818882298841197*^9}, {3.818882484197031*^9, 
   3.8188825534277687`*^9}, {3.818882598226076*^9, 3.818882640333383*^9}, {
   3.818888908281423*^9, 3.8188889278112392`*^9}, {3.818888967221086*^9, 
   3.8188889705271606`*^9}, 3.818937148818099*^9, {3.818937426136407*^9, 
   3.818937432570468*^9}, 3.818937475321644*^9, {3.8189381116498957`*^9, 
   3.8189381181436653`*^9}, {3.818938191935359*^9, 3.81893821437814*^9}, {
   3.818938430471423*^9, 3.81893848432169*^9}, {3.818942234753622*^9, 
   3.8189422427853727`*^9}, {3.818942330862174*^9, 3.8189423585351887`*^9}, {
   3.8189426798613234`*^9, 3.818942680455488*^9}, {3.819401928262059*^9, 
   3.819401933906714*^9}, 3.819402070631228*^9, {3.8196584166133757`*^9, 
   3.819658448313555*^9}, {3.819658480457584*^9, 3.819658481309371*^9}, {
   3.819990666528487*^9, 3.8199906927417717`*^9}, {3.8203237482832317`*^9, 
   3.8203237674806957`*^9}, {3.820358685996331*^9, 3.8203587119232883`*^9}, {
   3.820404758061131*^9, 3.82040477799226*^9}, {3.8204910736367817`*^9, 
   3.8204910868068*^9}, {3.820508610263049*^9, 3.820508632895162*^9}, {
   3.820525505875718*^9, 3.8205255064315767`*^9}, {3.820602072907857*^9, 
   3.820602085370233*^9}, {3.82061788546869*^9, 3.820617899424199*^9}, {
   3.8213084629722853`*^9, 3.821308473285358*^9}},
 CellLabel->
  "In[119]:=",ExpressionUUID->"a370cc62-5d4d-49d4-8ab4-8c404a66b998"],

Cell[BoxData["$Aborted"], "Output",
 CellChangeTimes->{3.820323829714547*^9, 3.820323902425563*^9, 
  3.82034451471565*^9, 3.82035866090021*^9, 3.8204047463009*^9, 
  3.8204910586604223`*^9, 3.8205086048730087`*^9, 3.820525494356079*^9, 
  3.820602050677581*^9, 3.820617857192227*^9, 3.820789098659465*^9, 
  3.821366831005777*^9},
 CellLabel->
  "Out[119]=",ExpressionUUID->"71840eeb-460b-4450-8987-1a9d7b230254"]
}, Open  ]]
},
WindowSize->{1115, 988},
WindowMargins->{{-7, Automatic}, {Automatic, 0}},
Magnification:>1.05 Inherited,
FrontEndVersion->"12.0 for Microsoft Windows (64-bit) (2019\:5e744\:67088\
\:65e5)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 8634, 225, 339, "Input",ExpressionUUID->"398a21c6-c9a3-4ee1-8c7d-444260a91e72"],
Cell[9195, 247, 248, 5, 49, "Input",ExpressionUUID->"71db4609-c829-4834-b1c4-1c1feed1af1c"],
Cell[9446, 254, 325, 7, 49, "Input",ExpressionUUID->"e1070aa4-d539-4b8c-b081-7402b6b16374"],
Cell[9774, 263, 427, 11, 49, "Input",ExpressionUUID->"d70282d0-7995-4523-b08c-3a45a4e5f7a1"],
Cell[10204, 276, 458, 12, 30, "Input",ExpressionUUID->"58628041-6b46-49e7-b651-4abb3f7a373a"],
Cell[10665, 290, 1853, 57, 75, "Input",ExpressionUUID->"e80a3ff6-3336-4d4f-a00d-95c9256a2780"],
Cell[12521, 349, 455, 12, 49, "Input",ExpressionUUID->"4c998990-0f46-443b-ac3a-a12e4bd34016"],
Cell[12979, 363, 3069, 96, 85, "Input",ExpressionUUID->"68b8b8fb-cf6f-4209-98ef-b208f8ed33c4"],
Cell[16051, 461, 1619, 49, 154, "Input",ExpressionUUID->"d2601902-6df0-4653-a0ef-ed3cfe9da0ae"],
Cell[17673, 512, 2971, 77, 239, "Input",ExpressionUUID->"ad538921-e09d-4dd7-a1d4-c57ba0f92e3f"],
Cell[20647, 591, 1502, 42, 157, "Input",ExpressionUUID->"19b11aa4-8a24-438f-b35a-d67d7ad30830"],
Cell[22152, 635, 1154, 34, 135, "Input",ExpressionUUID->"d0e2b5ef-98e0-4c6f-81b9-ad4ff82e8a0d"],
Cell[23309, 671, 3630, 97, 277, "Input",ExpressionUUID->"69850235-1da2-4e23-befc-4c70c9199bbb"],
Cell[26942, 770, 632, 14, 49, "Input",ExpressionUUID->"6c01a0ed-c53e-4a04-ac7b-64a4bf5a41d4"],
Cell[CellGroupData[{
Cell[27599, 788, 544, 11, 49, "Input",ExpressionUUID->"d447adda-a031-4fc8-8a5b-5fc66c8cc0b8"],
Cell[28146, 801, 1588, 38, 58, "Output",ExpressionUUID->"4174b2cb-305b-4658-bc33-7c5eec2aeeba"]
}, Open  ]],
Cell[29749, 842, 1068, 27, 90, "Input",ExpressionUUID->"4f3c4822-5dd2-4a93-9e50-e5e37aa00ae9"],
Cell[30820, 871, 2684, 53, 90, "Input",ExpressionUUID->"64b89cc3-7eae-45f0-9067-6ee016e76151"],
Cell[33507, 926, 1760, 50, 131, "Input",ExpressionUUID->"8be7cf6b-277d-4984-9bf5-23ee2193efa9"],
Cell[35270, 978, 269, 7, 49, "Input",ExpressionUUID->"890cd635-1b64-4a81-bc66-96214fa89eba"],
Cell[35542, 987, 1230, 24, 49, "Input",ExpressionUUID->"d96a3542-e0b4-4569-9ed6-fb27bc6953a4"],
Cell[36775, 1013, 8492, 203, 406, "Input",ExpressionUUID->"3f07b098-bae9-47b3-8b14-ad5b359acfb3"],
Cell[CellGroupData[{
Cell[45292, 1220, 292, 6, 30, "Input",ExpressionUUID->"5b85f204-51ae-41d9-bab1-f01ceae1173e"],
Cell[45587, 1228, 620, 11, 35, "Output",ExpressionUUID->"1775c23d-fba5-4c50-b85f-69933b890bf3"]
}, Open  ]],
Cell[46222, 1242, 953, 18, 49, "Input",ExpressionUUID->"3188ce0c-ef67-443c-986a-8db52dd3481a"],
Cell[47178, 1262, 1080, 22, 49, "Input",ExpressionUUID->"bdd7f63d-1c86-4ce1-a293-c5201e0c6a74"],
Cell[48261, 1286, 255, 6, 49, "Input",ExpressionUUID->"d2252a81-aac4-4e59-9d1c-3f7bec4d2273"],
Cell[CellGroupData[{
Cell[48541, 1296, 689, 14, 71, "Input",ExpressionUUID->"052d9d16-3cfe-4c04-b83b-ed3d9d4f54c0"],
Cell[49233, 1312, 539, 10, 35, "Output",ExpressionUUID->"91f4d2c4-7ed9-4e5c-afe9-ceec0a3e1658"],
Cell[49775, 1324, 16859, 278, 259, 171, 4, "CachedBoxData", "BoxData", "Output",ExpressionUUID->"4f6ad3a5-47c7-4c64-835e-8b8f959bf087"]
}, Open  ]],
Cell[CellGroupData[{
Cell[66671, 1607, 3160, 62, 217, "Input",ExpressionUUID->"a370cc62-5d4d-49d4-8ab4-8c404a66b998"],
Cell[69834, 1671, 414, 7, 35, "Output",ExpressionUUID->"71840eeb-460b-4450-8987-1a9d7b230254"]
}, Open  ]]
}
]
*)

